Vector Calculus 1 



The first rule in understanding vector calculus is draw lots of pictures. This subject can become rather 
abstract if you let it, but try to visualize all the manipulations. Try a lot of special cases and explore 
them. Keep relating the manipulations to the underlying pictures and don't get lost in the forest of 
infinite series. Along with the pictures, there are three types of derivatives, a couple of types of integrals, 
and some theorems relating them. 

9.1 Fluid Flow 

When water or any fluid moves through a pipe, what is the relationship between the motion of the 
fluid and the total rate of flow through the pipe (volume per time)? Take a rectangular pipe of sides a 
and b with fluid moving at constant speed through it and with the velocity of the fluid being the same 
throughout the pipe. It's a simple calculation: In time At the fluid moves a distance vAt down the 
pipe. The cross-section of the pipe has area A = ab, so the volume that move past a given flat surface 
is AV = AvAt. The flow rate is the volume per time, AV/At = Av. (The usual limit as At — > 
isn't needed here.) 




vAt ( a ) vAt (b) 



Just to make the problem look a little more involved, what happens to the result if I ask for the 
flow through a surface that is tilted at an angle to the velocity. Do the calculation the same way as 
before, but use the drawing (b) instead of (a). The fluid still moves a distance vAt, but the volume 
that moves past this flat but tilted surface is not its new (bigger) area A times vAt. The area of a 
parallelogram is not the product of its sides and the volume of a parallelepiped is not the area of a base 
times the length of another side. 




vAt 



The area of a parallelogram is the length of one side times the perpendicular distance from 
that side to its opposite side. Similarly the volume of a parallelepiped is the area of one side times 
the perpendicular distance from that side to the side opposite. The perpendicular distance is not the 
distance that the fluid moved (vAt). This perpendicular distance is smaller by a factor cos a, where 
a is the angle that the plane is tilted. It is most easily described by the angle that the normal to the 
plane makes with the direction of the fluid velocity. 

AV = Ah = A(vAt) cos a 

The flow rate is then AV/At = Ay cos a. Introduce the unit normal vector n, then this expression 
can be rewritten in terms of a dot product, 

— * 

Avcosa = Av-n = A-v (9.1) 
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where a is the angle between the direction of the fluid velocity and the normal to the area. 

This invites the definition of the area itself as a vector, and that's what I 
wrote in the final expression. The vector A is a notation for An, and defines the \ C x D 
area vector. If it looks a little odd to have an area be a vector, do you recall the 
geometric interpretation of a cross product? That's the vector perpendicular to 
two given vectors and it has a magnitude equal to the area of the parallelogram 
between the two vectors. It's the same thing. 

General Flow, Curved Surfaces 

The fluid velocity will not usually be a constant in space. It will be some function of position. The 
surface doesn't have to be flat; it can be cylindrical or spherical or something more complicated. How 
do you handle this? That's why integrals were invented. 

The idea behind an integral is that you will divide a complicated thing into small pieces and add 
the results of the small pieces to estimate the whole. Improve the estimation by making more, smaller 
pieces, and in the limit as the size of the pieces goes to zero get an exact answer. That's the procedure 
to use here. 

The concept of the surface integral is that you systematically divide a surface into a number (N) 
of pieces (k = 1, 2, . . .N). The pieces have area AA k and each piece has a unit normal vector n k . 
Within the middle of each of these areas the fluid has a velocity v k . This may not be a constant, but as 
usual with integrals, you pick a point somewhere in the little area and pick the v there; in the limit as 
all the pieces of area shrink to zero it won't matter exactly where you picked it. The flow rate through 
one of these pieces is Eq. (9.1), v k -n^AA^, and the corresponding estimate of the total flow through 
the surface is, using the notation AA k = n k AA k , 

N 

fc=i 

This limit as the size of each piece is shrunk to zero and correspondingly the number of pieces goes to 
infinity is the definition of the integral 

r - 

v-dA= lim ^v k -AA k (9.2) 
' k k=l 



Example of Flow Calculation 

In the rectangular pipe above, suppose that the flow exhibits shear, rising from zero at the bottom to 
Vo at the top. The velocity field is 

v(x,y,z) = v x (y)x = v^x (9.3) 

The origin is at the bottom of the pipe and the y-coordinate is measured upward from the origin. What 
is the flow rate through the area indicated, tilted at an angle <p from the vertical? The distance in and 
out of the plane of the picture (the 2-axis) is the length a. Can such a fluid flow really happen? Yes, 
real fluids such as water have viscosity, and if you construct a very wide pipe but not too high, and 
leave the top open to the wind, the horizontal wind will drag the fluid at the top with it (even if not as 
fast). The fluid at the bottom is kept at rest by the friction with the bottom surface. In between you 
get a gradual transition in the flow that is represented by Eq. (9.3). 
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Now to implement the calculation of the flow rate: 

Divide the area into N pieces of length Ai^ along the slant. 

The length in and out is a so the piece of area is AA^ = aA£^. 

The unit normal is % = xcos(f) — ysin0. (It happens to be independent of the index k, but that's 
special to this example.) 

The velocity vector at the position of this area is v = v^xy^/b. 

Put these together and you have the piece of the flow rate contributed by this area. 



Aflowfc = v ■ AAfc = v o 



b 

Vk 



x ■ a At], {x cos (J) — y sin 

£fc COS0 



Vo^aAi^ cos0 = t>o- 
o c 

In the last line I put all the variables in terms of £, using y = £cos< 
Now sum over all these areas and take a limit. 



lim 

A4-»o 



5> 



ifc COSl 



-aAlf- cos0 



bj cos <f> q 

divoj-icos 2 



aA£)~ cos (j) 



a? 2 2 j. 
V °b ~2 C ° S ^ 



1 COS (f> 



2b \cos0 



cos 2 <j) = v o 



ab 



This turns out to be independent of the orientation of the plane; the parameter is absent from the 
result. If you think of two planes, at angles <p\ and (f>2, what flows into one flows out of the other. 
Nothing is lost in between. 

Another Flow Calculation 

Take the same sort of fluid flow in a pipe, but make it a little more complicated. Instead of a flat 
surface, make it a cylinder. The axis of the cylinder is in and out in the picture and its radius is half 
the width of the pipe. Describe the coordinates on the surface by the angle 6 as measured from the 
midline. That means that — n/2 < 9 < n/2. Divide the surface into pieces that are rectangular strips 
of length a (in and out in the picture) and width bAO^/2. (The radius of the cylinder is b/2.) 

AAfr = a^-AOfr, and % = x cos Q\. + y sin 9j, (9.4) 



x 







k+l 




x 




<(b/2) sm6 



The velocity field is the same as before, v(x,y,z) = Voxy/b, so the contribution to the flow rate 
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through this piece of the surface is 



Vk ■ A A k = v yrr ■ a | A0 k n k 



The value of y k at the angle 6 k is 

b b . „ Vu 1 r ./in 

yfc = 2 + 2 sm fc ' so y = 2[ 1 + sin ^ 

Put the pieces together and you have 

v ^ [l + sin k ] x ■ a | A6» fc [x cos 6> fc + y sin fc ] = w ^ [l + sin 6> fc ] a | A# fc cos 6> fc 
The total flow is the sum of these over k and then the limit as A6 k — > 0. 

1 b r 71 "/ 2 1 6 

lim Vo-\l + sm9 k ]a-A9 k cos9 k = / Urj- [l + sin^la- d9 cos9 

Finally you can do the two terms of the integral: Look at the second term first. You can of course 
start grinding away and find the right trigonometric formula to do the integral, OR, you can sketch a 
graph of the integrand, sin#cos#, on the interval — 7r/2 < 6 < 7r/2 and write the answer down by 
inspection. The first part of the integral is 



ab W 2 . ab . . 
Vq— cos V = Vo — sin V 



4 J-k/2 ' 1 



tt/2 



ab 



-Tt/2 2 



And this is the same result as for the flat surface calculation. I set it up so that the two results are the 
same; it's easier to check that way. Gauss's theorem of vector calculus will guarantee that you get the 
same result for any surface spanning this pipe and for this particular velocity function. 

9.2 Vector Derivatives 

I want to show the underlying ideas of the vector derivatives, divergence and curl, and as the names 
themselves come from the study of fluid flow, that's where I'll start. You can describe the flow of a 
fluid, either gas or liquid or anything else, by specifying its velocity field, v(x,y,z) = v(f). 

/ / / - C 



For a single real-valued function of a real variable, it is often too complex to capture all the 
properties of a function at one glance, so it's going to be even harder here. One of the uses of ordinary 
calculus is to provide information about the local properties of a function without attacking the whole 
function at once. That is what derivatives do. If you know that the derivative of a function is positive 
at a point then you know that it is increasing there. This is such an ordinary use of calculus that you 
hardly give it a second thought (until you hit some advanced calculus and discover that some continuous 
functions don't even have derivatives). The geometric concept of derivative is the slope of the curve 
at a point — the tangent of the angle between the rr-axis and the straight line that best approximates 
the curve at that point. Going from this geometric idea to calculating the derivative takes some effort. 
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How can you do this for fluid flow? If you inject a small amount of dye into the fluid at some 
point it will spread into a volume that depends on how much you inject. As time goes on this region 
will move and distort and possibly become very complicated, too complicated to grasp in one picture. 




There must be a way to get a simpler picture. There is. Do it in the same spirit that you 
introduce the derivative, and concentrate on a little piece of the picture. Inject just a little bit of dye 
and wait only a little time. To make it explicit, assume that the initial volume of dye forms a sphere of 
(small) volume V and let the fluid move for a little time. 

1. In a small time At the center of the sphere will move. 

2. The sphere can expand or contract, changing its volume. 

3. The sphere can rotate. 

4. The sphere can distort. 

Div, Curl, Strain 

The first one, the motion of the center, tells you about the velocity at the center of the sphere. It's 
like knowing the value of a function at a point, and that tells you nothing about the behavior of the 
function in the neighborhood of the point. 

The second one, the volume, gives new information. You can simply take the time derivative 
dV/dt to see if the fluid is expanding or contracting; just check the sign and determine if it's positive 
or negative. But how big is it? That's not yet in a useful form because the size of this derivative 
will depend on how much the original volume is. If you put in twice as much dye, each part of the 
volume will change and there will be twice as much rate of change in the total volume. Divide the time 
derivative by the volume itself and this effect will cancel. Finally, to get the effect at one point take 
the limit as the volume approaches a point. This defines a kind of derivative of the velocity field called 
the divergence. 

lim — — j- = divergence of v (9.5) 
V-x) V dt 

This doesn't tell you how to compute it any more than saying that the derivative is the slope tells you 
how to compute an ordinary* derivative. I'll have to work that out. 

But first look at the third way that the sphere can move: rotation. Again, if you take a large 
object it will distort a lot and it will be hard to define a single rotation for it. Take a very small sphere 
instead. The time derivative of this rotation is its angular velocity, the vector uj. In the limit as the 
sphere approaches a point, this tells me about the rotation of the fluid in the immediate neighborhood 
of that point. If I place a tiny paddlewheel in the fluid, how will it rotate? 

2w = curlof{7 (9.6) 

The factor of 2 is for later convenience. 

After considering expansion and rotation, the final way that the sphere can change is that it 
can alter its shape. In a very small time interval, the sphere can slightly distort into an ellipsoid. This 
will lead to the mathematical concept of the strain. This is important in the subject of elasticity and 



* Can you start from the definition of the derivative as a slope, use it directly with no limits, and 
compute the derivative of x 2 with respect to x, getting 2x7 It can be done. 
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viscosity, but I'll put it aside for now save for one observation: how much information is needed to 
describe whatever it is? The sphere changes to an ellipsoid, and the first question is: what is the 
longest axis and how much stretch occurs along it — that's the three components of a vector. After 
that what is the shortest axis and how much contraction occurs along it? That's one more vector, 
but you need only two new components to define its direction because it's perpendicular to the long 
axis. After this there's nothing left. The direction of the third axis is determined and so is its length if 
you assume that the total volume hasn't changed. You can assume that is so because the question of 
volume change is already handled by the divergence; you don't need it here too. The total number of 
components needed for this object is 2 + 3 = 5. It comes under the heading of tensors. 

9.3 Computing the divergence 

Now how do you calculate these? I'll start with the simplest, the divergence, and compute the time 
derivative of a volume from the velocity field. To do this, go back to the definition of derivative: 

dt At^o At v ' 



"A .,--< 




vAt-~ : '"" 



Pick an arbitrary surface to start with and see how the volume changes as the fluid moves, 
carrying the surface with it. In time At a point on the surface will move by a distance vAt and it will 
carry with it a piece of neighboring area AA. This area sweeps out a volume. This piece of volume 
is not AA times vAt because the motion isn't likely to be perpendicular to the surface. It's only the 
component of the velocity normal to the surface that contributes to the volume swept out. Use n to 
denote the unit vector perpendicular to AA, then this volume is AAh-v At. This is the same as the 
calculation for fluid flow except that I'm interpreting the picture differently. 

If at a particular point on the surface the normal n is more or less in the direction of the velocity 
then this dot product is positive and the change in volume is positive. If it's opposite the velocity then 
the change is negative. The total change in volume of the whole initial volume is the sum over the 
entire surface of all these changes. Divide the surface into a lot of pieces AAi with accompanying unit 
normals n{, then 

A Vtotal = AA i ™* ' ^ At 

i 

Not really. I have to take a limit before this becomes an equality. The limit of this as all the AAi — > 
defines an integral 

AVtotai = <t> dAn-vAt 



and this integral notation is special; the circle through the integral designates an integral over the whole 
closed surface and the direction of n is always taken to be outward. Finally, divide by At and take the 
limit as At approaches zero. 

d ^ = jdAh-v (9.8) 

The v-hdA is the rate at which the area dA sweeps out volume as it is carried with the fluid. Note: 
There's nothing in this calculation saying that I have to take the limit as V — > 0; it's a perfectly 
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general expression for the rate of change of volume in a surface being carried with the fluid. It's also a 
completely general expression for the rate of flow of fluid through a fixed surface as the fluid moves past 
it. I'm interested in the first interpretation for now, but the second is just as valid in other contexts. 

Again, use the standard notation in which the area vector combines the unit normal and the area: 
dA = ndA. 

ldV .. 1 



divergence of v 



lim 

V^o V dt 



lim 

v^o V 



■dA 



(9.9) 



If the fluid is on average moving away from a point then the divergence there is positive. It's diverging. 
The Divergence as Derivatives 

This is still a long way from something that you can easily compute. I'll first go through a detailed 
analysis of how you turn this into a simple result, and I'll then go back to try to capture the essence 
of the derivation so you can see how it applies in a wide variety of coordinate systems. At that point 
I'll also show how to get to the result with a lot less algebra. You will see that a lot of the terms that 
appear in this first calculation will vanish in the end. It's important then to go back and see what was 
really essential to the calculation and what was not. As you go through this derivation then, try to 
anticipate which terms are going to be important and which terms are going to disappear. 

Express the velocity in rectangular components, v x x + v y y + v z z . For the small volume, choose 
a rectangular box with sides parallel to the axes. One corner is at point (xo,yo,Zo) and the opposite 
corner has coordinates that differ from these by (Ax, Ay, Az). Expand everything in a power series 
about the first corner as in section 2.5. Instead of writing out (xq, yo, Zo) every time, I'll abbreviate it 
by (o). 



Az 



x 



(x ,y ,z )Ax 



\^A 



v x (x,y,z) = 



Vx(o) + (x-X ) 

dv. 



a#(o) + G/-w>)-^(o) 



+ [x - x )(y - yo)^^(o) + • • 



(9.10) 



There are six integrals to do, one for each face of the box, and there are three functions, v x , v y , and 
v z to expand in three variables x, y, and z. Don't Panic. A lot of these are zero. If you look at the 
face on the right in the sketch you see that it is parallel to the y-z plane and has normal n = x. When 
you evaluate v-n, only the v x term survives; flow parallel to the surface (v y , v z ) contributes nothing 
to volume change along this part of the surface. Already then, many terms have simply gone away. 

Write the two integrals over the two surfaces parallel to the y-z plane, one at Xq and one at 
x + Ax. 




ryo+Ay rzo+Az ryo+Ay rz +Az 

I dy dzv x (x + Ax,y,z) - / dy dzv x (x ,y,z) 

yo Jz Jyo Jz 



The minus sign comes from the dot product because n points left on the left side. You can evaluate 
these integrals by using their power series representations, and though you may have an infinite number 
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of terms to integrate, at least they're all easy. Take the first of them: 

«y Q +ky fZ +Az 



\ dy dz 

Jy Jzq 



V x {o) + 



= v x ( )AyAz + (Ax)-^( )AyAz+ 

i(A,) 2 A^( ) + '-(AzfAy^io) + i(Ax) 2 g|(o)A y A, + • • • 

Now look at the second integral, the one that you have to subtract from this one. Before plunging in to 
the calculation, stop and look around. What will cancel; what will contribute; what will not contribute? 
The only difference is that this is now evaluated at Xq instead of at Xq + Ax. The terms that have 
Ax in them simply won't appear this time. All the rest are exactly the same as before. That means 
that all the terms in the above expression that do not have a Ax in them will be canceled when you 
subtract the second integral. All the terms that do have a Ax will be untouched. The combination of 
the two integrals is then 

(A*)^( )A^ + l -(Axf^( )AyAz + I(Ax)^( )(Ay) 2 A^ + • • • 

Two down four to go, but not really. The other integrals are the same except that x becomes y 
and y becomes z and z becomes x. The integral over the two faces with y constant are then 

(jI) i CJ 11 

(Ay)^( )AzAx + -(Ayf^-( )AzAx + ■■■ 

and a similar expression for the final two faces. The definition of Eq. (9.9) says to add all three of these 
expressions, divide by the volume, and take the limit as the volume goes to zero. V = AxAyAz, and 
you see that this is a common factor in all of the terms above. Cancel what you can and you have 

dv x dv y dv x 1 ,d 2 v x . l tA ,d 2 v y 1 ,d 2 v z 

7^ (o) + W (o) + ~dx- [o) + 2 (Ax) a^ (o) + 2 {Ay) ^ io) + 2 (A ^ (o) + • • • 

In the limit that the all the Ax, Ay, and Az shrink to zero the terms with a second derivative vanish, as 
do all the other higher order terms. You are left then with a rather simple expression for the divergence. 

,. -, dv x dv v dv x 

divergence of v = divt> = h -ft- + -ft- (9.11) 

ox ay ox 

This is abbreviated by using the differential operator V, "del." 

~ - 9 „ d „ d i , 

V = X d-x +y dy +Z dz- (9 ' 12) 

Then you can write the preceding equation as 

divergence of v = divv = V ■ v (9.13) 



9 — Vector Calculus 1 



9 



The symbol V will take other forms in other coordinate systems. 

Now that you've waded through this rather technical set of manipulations, is there an easier way? 
Yes but, without having gone through the preceding algebra you won't be able to see and to understand 
which terms are important and which terms are going to cancel or otherwise disappear. When you need 
to apply these ideas to something besides rectangular coordinates you have to know what to keep and 
what to ignore. Once you know this, you can go straight to the end of the calculation and write down 
those terms that you know are going to survive, dropping the others. This takes practice. 

Simplifying the derivation 

In the long derivation of the divergence, the essence is that you find v ■ n on one side of the box (maybe 
take it in the center of the face), and multiply it by the area of that side. Do this on the other side, 
remembering that n isn't in the same direction there, and combine the results. Do this for each side 
and divide by the volume of the box. 

[v x (x + Ax,y + Ay/2, z + Az/2)AyAz 

-v x (x , y + Ay/2, z + Az/2)AyAz] 4- (AxAyAz) (9.14) 

the Ay and Az factors cancel, and what's left is, in the limit Ax — > 0, the derivative dv x /dx. 

I was careful to evaluate the values of v x in the center of the sides, but you see that it didn't 
matter. In the limit as all the sides go to zero I could just as easily taken the coordinates at one corner 
and simplified the steps still more. Do this for the other sides, add, and you get the result. It all looks 
very simple when you do it this way, but what if you need to do it in cylindrical coordinates? 

-rA<fi 

■^r^^-^Az 

Ar 

When everything is small, the volume is close to a rectangular box, so its volume is V = 
(Ar)(Az)(rA(f)). Go through the simple version for the calculation of the surface integral. The top 
and bottom present nothing significantly different from the rectangular case. 

dv 

[v z (r , 0o, z + Az) - v z (r , 0o, z )] (Ar)(r o A0) -=- r Ar A(f)Az — ► 

The curved faces of constant r are a bit different, because the areas of the two opposing faces 
aren't the same. 

1 (j(t"v ) 

[v r (r + Ar,(f) ,z )(r + Ar)A(f)Az -v r (r ,<f)o,Zo)roA(f)Az] ~r ArA(f)Az — ► - 

A bit more complex than the rectangular case, but not too bad. 

Now for the constant sides. Here the areas of the two faces are the same, so even though they 
are not precisely parallel to each other this doesn't cause any difficulties. 




[v<p(r , 0o + A(f), zq) - v z (r , O , z )] (Ar)(Az) ~ r Ar Acf)Az 



l dv 4 
r (90 



The sum of all these terms is the divergence expressed in cylindrical coordinates. 

^ ld(rv r ) .ldvj, dv z 

divv = -^- — - -| -f + ^— (9.15) 

r or r o<p oz 
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The corresponding expression in spherical coordinates is found in exactly the same way, prob- 
lem 9.4. 

,. _ 1 d(r 2 v r ) 1 d(sm9v s ) 1 dru 
r A or rsmO o9 rsmv o<p 

These are the three commonly occurring coordinates system, though the same simplified method 
will work in any other orthogonal coordinate system. The coordinate system is orthogonal if the surfaces 
made by setting the value of the respective coordinates to a constant intersect at right angles. In the 
spherical example this means that a surface of constant r is a sphere. A surface of constant 9 is a 
half-plane starting from the 2-axis. These intersect perpendicular to each other. If you set the third 
coordinate, 0, to a constant you have a cone that intersects the other two at right angles. Look back 
to section 8.8. 

9.4 Integral Representation of Curl 

The calculation of the divergence was facilitated by the fact that the equation (9.5) could be manipulated 
into the form of an integral, Eq. (9.9). Is there a similar expression for the curl? Yes. 

curl-u = lim — <p dA x v (9.17) 
V^o V J 

For the divergence there was a logical and orderly development to derive Eq. (9.9) from (9.5). Is there 
a similar intuitively clear path here? I don't know of one. The best that I can do is to show that it 
gives the right answer. 

And what's that surface integral doing with a x instead of a ■? No mistake. Just replace the 
dot product by a cross product in the definition of the integral. This time however you have to watch 
the order of the factors. 

,w 

ndA 




To verify that this does give the correct answer, use a vector field that represents pure rigid body 
rotation. You're going to take the limit as AV — > 0, so it may as well be uniform. The velocity field 
for this is the same as from problem 7.5. 

v = oj x r (9.18) 

To evaluate the integral use a sphere of radius R centered at the origin, making n = f . You also need 
the identity A x (B x C) = B(A-C) - C(A-B). 

dAx (u x f) =u(dA-f) -f(u-dA) (9.19) 
Choose a spherical coordinate system with the z-axis along uj. 



— * —> 

dA = ndA = fdA, and Q ■ dA = u dA cos 9 
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j> dA xv = j> QR dA — ruj dA cos 9 



uR4irR 2 -u / R 2 sm9d9 / d<j> z R cos 9 cos 9 
Jo Jo 

uR AnR 2 - uz 2itR 3 / sin 9d9 cos 2 9 



o 

/ 1 2 
cos 2 9dcos9 = ooR AttR 2 - loz 2ttR 3 ■ - 
i 3 

= lu-ttR 3 
3 

Divide by the volume of the sphere and you have 2u as promised. In the first term on the first line of 
the calculation, tuR is a constant over the surface so you can pull it out of the integral. In the second 
term, r has components in the x, y, and z directions; the first two of these integrate to zero because 
for every vector with a positive ^-component there is one that has a negative component. Same for y. 
All that is left of r*is 2Rcos9. 

The Curl in Components 

With the integral representation, Eq. (9.17), available for the curl, the process is much like that for 
computing the divergence. Start with rectangular of course. Use the same equation, Eq. (9.10) and 
the same picture that accompanied that equation. With the experience gained from computing the 
divergence however, you don't have to go through all the complications of the first calculation. Use the 
simpler form that followed. 

In Eq. (9.14) you have v ■ AA = v x AyAz on the right face and on the left face. This time 
replace the dot with a cross (in the right order). 

On the right, 

AA xv = AyAz x x v(x + Ax, y + Ay/2, z + Az/2) (9.20) 

On the left it is 

AA xv = AyAz x x v(x , y + Ay/2, z + Az/2) (9.21) 

When you subtract the second from the first and divide by the volume, AxAyAz, what is left is (in 
the limit Ax — > 0) a derivative. 

- v{xp + Ax, 7/0, go) - v(x , 2/q, z ) . dv 
Ax ox 

„ f.dv x „dv y „dv z 
= x x x^t — h y-^r- + z-ft- 
\ ox ox ox 

__ *dv y ^dv z 

Z dx y dx 

Similar calculations for the other four faces of the box give results that you can get simply by 
changing the labels: x— > y — > z — > x, a cyclic permutation of the indices. The result can be expressed 
most succinctly in terms of V. 

curlt> = Vxtf (9.22) 

In the process of this calculation the normal vector x was parallel on the opposite faces of the 
box (except for a reversal of direction). Watch out in other coordinate systems and you'll see that this 
isn't always true. Just draw the picture in cylindrical coordinates and this will be clear. 



9 — Vector Calculus 1 



12 



9.5 The Gradient 

The gradient is the closest thing to an ordinary derivative here, taking a scalar-valued function into a 
vector field. The simplest geometric definition is "the derivative of a function with respect to distance 
along the direction in which the function changes most rapidly," and the direction of the gradient vector 
is along that most-rapidly-changing direction. If you're dealing with one dimension, ordinary real-valued 
functions of real variables, the gradient is the ordinary derivative. Section 8.5 has some discussion and 
examples of this, including its use in various coordinate systems. It is most conveniently expressed in 
terms of V. 

grad/ = V/ (9.23) 

The equations (8.15), (8.27), and (8.28) show the gradient (and correspondingly V) in three 
coordinate systems. 

„ - d „ d „ d 
rectangular: V = x^ — hy^ — h^Tr 
ox dy dz 

d -1 d d 
cylindrical: V = f^- + + z-*- (9.24) 

dr r d<p dz 

, . , „ „d ,19 7 1 d 
spherical: V = — \-"~^ + 



dr r d6 r sin 6 d<p 

In all nine of these components, the denominator (e.g. rsinOdcj)) is the element of displacement along 
the direction indicated. 

9.6 Shorter Cut for div and curl 

There is another way to compute the divergence and curl in cylindrical and rectangular coordinates. A 
direct application of Eqs. (9.13), (9.22), and (9.24) gets the result quickly. The major caution is that 
you have to be careful that the unit vectors are inside the derivative, so you have to differentiate them 
too. 

V ■ v is the divergence of v, and in cylindrical coordinates 

V '' 7= {^o\ + ^\§4> +2 §^) -irVr + facP + ZVz) (9.25) 

The unit vectors f, </>, and z don't change as you alter r or z. They do change as you alter 0. (except 
for z). 

df _ d(p _ dz _ dr _ d(p _ dz _ dz _ 
dr dr dr dz dz dz dcj) 
Next come df/d<p and d(p/d(j). This is problem 8.20. You can do this by first showing that 

f = x cos (f) + y sin <j> and <f> = — xsmcf) + y cos0 (9.27) 

and differentiating with respect to (p. This gives 

df/d<f) = 4>, and d$/d<p = -r (9.28) 

Put these together and you have 

„ _ dv r t 1 d 7 x dv z 

V - V = W + ^rd6 ( - rVr + ^ ) + ^z- 



dr r \ r dcf) dcj) ) dz 
dv r 1 1 dvs dv z 

-x L + -Vr + -^x + ^- 9.29 
dr r r dm dz 
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This agrees with equation (9.15). 

Similarly you can use the results of problem 8.15 to find the derivatives of the corresponding 
vectors in spherical coordinates. The non-zero values are 

4t = 4> sin 9 -tt = $ cos 6 -^r = —f sin # — cos 6 

d(p d(p d(p 

dr % dd , 

m =B dS = - r (9 - 30) 

The result is for spherical coordinates 

V .tf = + 1 d(smOv e ) + _J_dv l 

r 2 dr rsinO d6 rsinO d<p 

The expressions for the curl are, cylindrical: 

V x v = f t 1 - 9 ^ - ^ + 3>(*± - 9 M + 2 f 1 - 9 ^- - 1 - 9 ^] (9.32) 
\r d(f) dz J \ dz dr J \r dr r d<p J 

and spherical: 

_ ^ „ 1 fdismOvs) dv 6 
V x v = r — ' 

r s 

^^rsin^ dcj) r dr ) ^ r ^ dr d9 ) (9.33) 
9.7 Identities for Vector Operators 

Some of the common identities can be proved simply by computing them in rectangular components. 
These are vectors, and if you show that one vector equals another vector it doesn't matter that you 
used a simple coordinate system to demonstrate the fact. Of course there are some people who quite 
properly complain about the inelegance of such a procedure. They're called mathematicians. 

V-Vxv = VxV/ = V x V x v = V(V-v) - (V ■ V)v (9.34) 

There are many other identities, but these are the big three. 

v-dA= d 3 rV-v <k v-df= / V xv-dA (9.35) 



rsin# V d9 



1 dv r l<9(r^)\ -l fd{rv e ) dv r 



are the two fundamental integral relationships, going under the names of Gauss and Stokes. See chapter 
13 for the proofs of these integral relations. 

9.8 Applications to Gravity 

The basic equations to describe the gravitational field in Newton's theory are 



V ■ £ = -4vrGp, and Vxj = (9.36) 

In these equations, the vector field g is defined by placing a (very small) test mass m at a point and 
measuring the gravitational force on it. This force is proportional to m itself, and the proportionality 
factor is called the gravitational field g. The other symbol used here is p, and that is the volume mass 
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density, dm/dV of the matter that is generating the gravitational field. G is Newton's gravitational 
constant: G = 6.67 x 10 -11 N-m 2 /kg 2 . 




For the first example of solutions to these equations, take the case of a spherically symmetric 
mass that is the source of a gravitational field. Assume also that its density is constant inside; the 
total mass is M and it occupies a sphere of radius R. Whatever g is, it has only a radial component, 
g = g r f. Proof: Suppose it has a sideways component at some point. Rotate the whole system by 
180° about an axis that passes through this point and through the center of the sphere. The system 
doesn't change because of this, but the sideways component of g would reverse. That can't happen. 

The component g r can't depend on either 9 or (ft because the source doesn't change if you rotate 
it about any axis; it's spherically symmetric. 

g = 9r(r)f (9.37) 
Now compute the divergence and the curl of this field. Use Eqs. (9.16) and (9.33) to get 

1 d(r 2 g r \ 

V ■ q r (r)f = — — - and V x o r (r)f = 
r 2 dr 

The first equation says that the divergence of g is proportional to p. 

\&l = -^Gp (9.38) 
r 1 dr 

Outside the surface r = R, the mass density is zero, so this is 

1 d(r 2 g r ) . , . 2 ^ , C 

^ dr = °> implying rg r = C, and g r = ^ 

where C is some as yet undetermined constant. Now do this inside. 

1 d(r 2 g r ) 



r 2 dr 



-AnGp , where p = 3M/4itR 3 



This is 



d[rl ' lr) - -47cGp r 2 , so r 2 g r = -^Gp^ + C\ 



or g r (r) = --irGp r + 



4 C 
-7iGp r + — 
3 r z 

There are two constants that you have to evaluate: C and C. The latter has to be zero, because 
C'/r 2 — > 00 as r — > 0, and there's nothing in the mass distribution that will cause this. As for the 
other, note that g r must be continuous at the surface of the mass. If it isn't, then when you try to 
differentiate it in Eq. (9.38) you'll be differentiating a step function and you get an infinite derivative 
there (and the mass density isn't infinite there). 

4 C 
g r (R-) = --irGpo R = g r {R+) = 
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Solve for C and you have 



C = --nG Po R 3 



4 3M D3 
3 47r_R 3 



-CM 



Put this all together and express the density po m terms of M and i? to get 

R 

9r{r) = 



9i 




GM/r 2 (r > i?) 
GMr I R? (r < i?) 



(9.39) 



This says that outside the spherical mass distribution you can't tell what its radius R is. It 
creates the same gravitational field as a point mass. Inside the uniform sphere, the field drops to zero 
linearly toward the center. For a slight variation on how to do this calculation see problem 9.14. 

Non-uniform density 

The density of the Earth is not uniform; it's bigger in the center. The gravitational field even increases 
as you go down below the Earth's surface. What does this tell you about the density function? V ■ g = 
—AnGp remains true, and I'll continue to make the approximation of spherical symmetry, so this is 



1 d(r 2 g r ) _ dg r 2 
r 2 dr dr r^ r 71 p r 

That gravity increases with depth (for a little while) says 



(9.40) 



dgr 
dr 



-4iiGp(r) - -g r > 



Why > 0? Remember: g r is itself negative and r is measured outward. Sort out the signs. I can solve 
for the density to get 

1 

p(r) < 



2nGr 



9r 



At the surface, g r (R) = -GM/R 2 , so this is 



p(R)< 



M 2 3M 



2vri?3 " 3 " 4ttR3 ~ 3 Paverage 



The mean density of the Earth is 5.5gram/cm 3 , so this bound is 3.7gram/ctn 3 . Pick up a random 
rock. What is its density? 



9.9 Gravitational Potential 

The gravitational potential is that function V for which 

g = ~VV 



(9.41) 



That such a function even exists is not instantly obvious, but it is a consequence of the second of the 
two defining equations (9.36). If you grant that, then you can get an immediate equation for V by 
substituting it into the first of (9.36). 



V ■ g = -V ■ VV = -ArrGp, or V 2 V = 4itGp 



(9.42) 
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This is a scalar equation instead of a vector equation, so it will often be easier to handle. Apply it to 
the same example as above, the uniform spherical mass. 

The Laplacian, V 2 is the divergence of the gradient, so to express it in spherical coordinates, 
combine Eqs. (9.24) and (9.31). 

1 d f 2 dV\ , 1 d ( . n dV\ . 1 d 2 V 



V V ~ r 2 dr [ r dr ) + r 2 gin QQ [ S ' m 9 QQ ) + r 2 gin 2 Q d(j) 2 ( 9 " 43 ) 

Because the mass is spherical it doesn't change no matter how you rotate it so the same thing 
holds for the solution, V{r). Use this spherical coordinate representation of V 2 and for this case the 6 
and (f) derivatives vanish. 

I changed from d to d because there's now one independent variable, not several. Just as with Eq. (9.38) 
I'll divide this into two cases, inside and outside. 

n . , I d f 2 dV\ 2 dV _ 

Outside: — —{r 2 -^)=0, so r 2 -r- = C 



r 2 dr \ dr j dr 
Continue solving this and you have 

dV C C 

w = F2 ^V(r) = - T + D (r>R) (9.45) 

, . , 1 d / 2 dV\ t _ 2 dV A „ r 3 . 

Inside: ^ I r \ = 47rGp so r z = AirGpo— + C 

Continue, dividing by r 2 and integrating, 

V(r) = 4irGp j - — + D' (r < R) (9.46) 

There are now four arbitrary constants to examine. Start with C . It's the coefficient of l/r in the 
domain where r < R. That means that it blows up as r — > 0, but there's nothing at the origin to 
cause this. C = 0. Notice that the same argument does not eliminate C because (9.45) applies only 
for r > R. 

Boundary Conditions 

Now for the boundary conditions at r = R. There are a couple of ways to determine this. I find 
the simplest and the most general approach is to recognize that the equations (9.42) and (9.44) must 
be satisfied everywhere. That means not just outside, not just inside, but at the surface too. The 
consequence of this statement is the result* 

V is continuous at r = R dV/dr is continuous at r = R (9.47) 

Where do these continuity conditions come from? Assume for a moment that the first one is false, that 
V is discontinuous at r = R, and look at the proposition graphically. If V changes value in a very 
small interval the graphs of V, of dV/dr, and of d 2 V/dr 2 look like 



* Watch out for the similar looking equations that appear in electrostatics. Only the first of these 
equations holds there; the second must be modified by a dielectric constant. 
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V 



dV/dr 



d 2 V/dr 2 



The second derivative on the left side of Eq. (9.44) has a double spike that does not appear on 
the right side. It can't be there, so my assumption that V is discontinuous is false and V must be 
continuous. 

Assume next that V is continuous but its derivative is not. The graphs of V, of dV/dr, and of 
d 2 V/dr 2 then look like 




dV/dr 



d 2 V/dr 2 



The second derivative on the left side of Eq. (9.44) still has a spike in it and there is no such 
spike in the p on the right side. This is impossible, so dV/dr too must be continuous. 

Back to the Problem 

Of the four constants that appear in Eqs. (9.45) and (9.46), one is already known, C . For the rest, 



V(R-) = V(R+) 



is 



o 



c 

R 



+ D 



8nGpo 



R 



C_ 
R? 



These two equations determine two of the constants 

R 3 



C = AnGpc 



then 



D-D' = AnGp ^-+^Gp ^- 
b 6 



2TcGp R 2 



Put this together and you have 



V 



R 



V(r) 



inGp r 2 - 27rGp R 2 + D (r < R) 



T Ti 



Gp R 3 /r + D 



(r > R) 



(9.48) 



Did I say that the use of potentials is supposed to simplify the problems? Yes, but only the harder 
problems. The negative gradient of Eq. (9.48) should be g. Is it? The constant D can't be determined 
and is arbitrary. You may choose it to be zero. 
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Magnetic Boundary Conditions 

The equations for (time independent) magnetic fields are 

V x B = hqJ and V ■ B = (9.49) 

The vector J is the current density, the current per area, defined so that across a tiny area dA the 
current that flows through the area is dl = J ■ dA. (This is precisely parallel to Eq. (9.1) for fluid 
flow rate.) In a wire of radius R, carrying a uniform current /, the magnitude of J is I /irR 2 . These 
equations are sort of the reverse of Eq. (9.36). 

If there is a discontinuity in the current density at a surface such as the edge of a wire, will 
there be some sort of corresponding discontinuity in the magnetic field? Use the same type of analysis 
as followed Eq. (9.47) for the possible discontinuities in the potential function. Take the surface of 
discontinuity of the current density to be the x-y plane, z = and write the divergence equation 

dB x dB y dB z 




-I = o 

dx dy dz 

If there is a discontinuity, it will be in the z variable. Perhaps B x or B z is discontinuous at the x-y plane. 
The divergence equation has a derivative with respect to z only on B z . If one of the other components 
changes abruptly at the surface, this equation causes no problem — nothing special happens in the x 
or y direction. If B z changes at the surface then the derivative dB z /dz has a spike. Nothing else in 
the equation has a spike, so there's no way that you can satisfy the equation. Conclusion: The normal 
component of B is continuous at the surface. 
What does the curl say? 

X \dy dx) +y \dz dy) +Z \dx dy )-^{xJ x + yJ y + zJ z ) 

Derivatives with respect to x or y don't introduce a problem at the surface, all the action is again 
along z. Only the terms with a d /dz will raise a question. If B x is discontinuous at the surface, then 
its derivative with respect to z will have a spike in the y direction that has no other term to balance 
it. (J y has a step here but not a spike.) Similarly for B y . This can't happen, so the conclusion: The 

tangential component of B is continuous at the surface. 

What if the surface isn't a plane. Maybe it is a cylinder or a sphere. In a small enough region, 
both of these look like planes. That's why there is still* a Flat Earth Society. 

9.10 Index Notation 

In section 7.11 I introduced the summation convention for repeated indices. I'm now going to go over 
it again and emphasize its utility in practical calculations. 

When you want to work in a rectangular coordinate system, with basis vectors x, y, and z, it is 
convenient to use a more orderly notation for the basis vectors instead of just a sequence of letters of 
the alphabet. Instead, call them e\, §2, and £3. (More indices if you have more dimensions.) I'll keep 
the assumption that these are orthogonal unit vectors so that 

ei -e 2 = 0, e 3 -e 3 = 1, etc. 

More generally, write this in the compact notation 

8i-2i = fy = {j (i^-j ( 9 - 5 °) 



en.wikipedia.org/wiki/Flat_Earth_Society 
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The Kronecker delta is either one or zero depending on whether i = j or i ^ j, and this equation sums 
up the properties of the basis in a single, compact notation. You can now write a vector in this basis as 

— * 

A = Ai ei + A 2 e 2 + A 3 e 3 = A h 

The last expression uses the summation convention that a repeated index is summed over its range. 
When an index is repeated in a term, you will invariably have exactly two instances of the index; if you 
have three it's a mistake. 

When you add or subtract vectors, the index notation is 

A + B = C = A i e i + B i e i = (A i + B i )e i = C i e i or A i + B i = C i 

Does it make sense to write A + B = A^e^ + B k e k 7 Yes, but it's sort of pointless and confusing. You 
can change any summed index to any label you find convenient — they're just dummy variables. 

Aiii = A e e e = A m e m = A 1 e 1 + A 2 e 2 + A 3 e 3 

You can sometime use this freedom to help do manipulations, but in the example Ai&i + B k e k it's no 
help at all. 

Combinations such as 

E i + Fi or E k F k G % = H % or M ke D e = F k 

are valid. The last is simply Eq. (7.8) for a matrix times a column matrix. 

Ai + Bj = C k or E m F m G m or C k = AijBj 

have no meaning. 

You can manipulate the indices for your convenience as long as you follow the rules. 

Ai = BijCj is the same as A k = B kn C n or A^ = Bg^Cp 
The scalar product has a simple form in index notation: 

A-B = Aiti-Bj ej = AiBj e { ■ ij = AiBjSij = A { Bi (9.51) 

The final equation comes by doing one of the two sums (say j), and only the term with j = i survives, 
producing the final expression. The result shows that the sum over the repeated index is the scalar 
product in disguise. 

Just as the dot product of the basis vectors is the delta-symbol, the cross product provides 
another important index function, the alternating symbol. 

(9.52) 





6j ■ Cj 


x e k = 






<§2 x e 3 = 


yxz = e ± , 


so 


Cl2 3 


= <°i ■ e 2 x e 3 : 


e 2 x e 3 


= -e 3 x e 2 , 


so 


ei 3 2 


= — Cl2 3 = — 1 


e 3 ■ §i x e 2 = 


: e 3 -e 3 = 1, 


so 


e 3 i2 


= ei2 3 = 1 


ei ■ e 3 x e 3 


= ei ■ = 


so 


^i 33 


= 
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If the indices are a cyclic permutation of 123, (231 or 312), the alternating symbol is 1. 

If the indices are an odd permutation of 123, (132 or 321 or 213), the symbol is — 1. 

If any two of the indices are equal the alternating symbol is zero, and that finishes all the cases. The 

last property is easy to see because if you interchange any two indices the sign changes. If the indices 

are the same, the sign can't change so it must be zero. 

Use the alternating symbol to write the cross product itself. 

A x B = Ai§i x Bj&j and the /c-component is 

&k ' ^ ^ ~~ ^k ' ^-i-^j ^ ~~ ^kij-^-i-^j 

You can use the summation convention to advantage in calculus too. The V vector operator has 
components 

Vj or some people prefer di 
For unity of notation, use X\ = x and X2 = y and 23 = z. In this language, 

f3i = Vi is ^- = w- (9-53) 

OX OXi 

Note: This notation applies to rectangular component calculations only! (e^-ej = 5{j.) The generalization 
to curved coordinate systems will wait until chapter 12. 

,. ^ „ - r, dvi dv 2 dv 3 

divv = V-v = d i v i = 1 r A + 1 r A + 1 r A 9.54 

OX 1 OX2 OX3 

You should verify that d{Xj = 5ij. 

Similarly the curl is expressed using the alternating symbol. 

cuvlv = V x v becomes ^ijk^j v k = (curl-y^ (9.55) 

the i th components of the curl. 

An example of manipulating these object: Prove that curlgrad0 = 0. 

curlgrad^ = V x V(f> — v e^djd^ (9.56) 

You can interchange the order of the differentiation as always, and the trick here is to relabel the indices 
— that is a standard technique in this business. 

e ijkdjdk4> = eijkdkdj(fi = £ikjdjd k (p 

The first equation interchanges the order of differentiation. In the second equation, I call "j" "k" and 
I call "fc" "j". These are dummy indices, summed over, so this can't affect the result, but now this 
expression looks like that of Eq. (9.56) except that two (dummy) indices are reversed. The e symbol 
is antisymmetric under interchange of any of its indices, so the final expression is the negative of the 
expression in Eq. (9.56). The only number equal to minus itself is zero, so the identity is proved. 
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9.11 More Complicated Potentials 

The gravitational field from a point mass is g 
(p = —Gm/r. This satisfies 



-Gmf/r 2 , so the potential for this point mass is 



g = -V<P = -V- 



-Gm 



, d Gm 

dr r 



Gmf 



For several point masses, the gravitational field is the vector sum of the contributions from each mass. 
In the same way the gravitational potential is the (scalar) sum of the potentials contributed by each 
mass. This is almost always easier to calculate than the vector sum. If the distribution is continuous, 
you have an integral. 



0. 



total 



Gm^ 



or 



/ 



G dm 




This sort of very abbreviated notation for the sums and integrals is normal once you have done a lot of 
them, but when you're just getting started it is useful to go back and forth between this terse notation 
and a more verbose form. Expand the notation and you have 



/dun 
— ^7 
r _ r i 



(9.57) 



This is still not very explicit, so expand it some more. Let 



then 



f' = xx' + yy' + zz' and 
x,y,z) = —G / dx'dy'dz' p(x' ,y' , z') 



r = xx + yy + zz 
1 



y/[x - x'f + {y- y') 2 + {z- z'Y 



where p is the volume mass density so that dm = pdV = pd 3 r', and the limits of integration are such 
that this extends over the whole volume of the mass that is the source of the potential. The primed 
coordinates represent the positions of the masses, and the non-primed ones are the position of the point 
where you are evaluating the potential, the field point. The combination d ? 'r is a common notation for 
a volume element in three dimensions. 

For a simple example, what is the gravitational potential from a uniform thin rod? Place its 
center at the origin and its length = 2L along the z-axis. The potential is 



Gdm 



-G 



Xdz' 



^ X 2 +y 2 + ( z _ z ,y 



where A = M/2L is its linear mass density. This is an elementary integral. Let u = z' — z, and 

a = \/ x 2 + y 2 . 



= -GX 



L-z 



L-z 



du 
\J a 2 + u 2 



-GX / d9 = ~GX9 



u=L—z 



-L-z 
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where u = asinh#. Put this back into the original variables and you have 



-GX 



sinh 



-i 



L 



\Jx 2 + y 2 



+ sinh 




(9.58) 



The inverse hyperbolic function is a logarithm as in Eq. (1.4), so this can be rearranged and the terms 
combined into the logarithm of a function of x, y, and z, but the sinh - s are easier to work with so 
there's not much point. This is not too complicated a result, and it is far easier to handle than the 
vector field you get if you take its gradient. It's still necessary to analyze it in order to understand it 
and to check for errors. See problem 9.48. 



Exercises 

1 Prove that the geometric interpretation of a cross product is as an area. 

2 Start from a picture of C = A — B and use the definition and properties of the dot product to derive 
the law of cosines. (If this takes you more than about three lines, start over, and no components, just 
vectors.) 

3 Start from a picture of C = A — B and use the definition and properties of the cross product to 
derive the law of sines. (If this takes you more than a few lines, start over.) 

4 Show that A-BxC = AxB-C. Do this by drawing pictures of the three vectors and find the 
geometric meaning of each side of the equation, showing that they are the same (including sign). 

5 (a) If the dot product of a given vector F with every vector results in zero, show that F = 0. 
(b) Same for the cross product. 

6 From the definition of the dot product, and in two dimensions, draw pictures to interpret A -(B + C) 
and from there prove the distributive law: A • = A - B + A - C '. (Draw B + C tip-to-tail.) 

7 For a sphere, from the definition of the integral, what is § dAl What is § dAl 

8 What is the divergence of xxy + yyz + zzxl 

9 What is the divergence of fr sin 6 + 6r$ sin 6 cos + 0r cos 0? (spherical) 

10 What is the divergence of fr sin <p + <pz sin (p + zzrl (cylindrical) 

11 In cylindrical coordinates draw a picture of the vector field V = 0r 2 (2+cos 0) (for z = 0). Compute 
the divergence of V and indicate in a second sketch what it looks like. Draw field lines for v. 

12 What is the curl of the vector field in the preceding exercise (and indicate in a sketch what it's like). 

13 Calculate V-f. (a) in cylindrical and (b) in spherical coordinates. 

14 Compute diXj. 

15 Compute div curl v using index notation. 

16 Show that = |(z — — k)(k — i). 

17 Use index notation to derive V -(fv) = (V/) ■ V + /V ■ V. 

18 Write Eq. (8.6) in index notation, where (x,y) — > x^ and (x',y') — > x'-. 
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Problems 

9.1 Use the same geometry as that following Eq. (9.3), and take the velocity function to be V = 
xvoxy/b 2 . Take the bottom edge of the plane to be at (x,y) = (0,0) and calculate the flow rate. 
Should the result be independent of the angle 0? Sketch the flow to understand this point. Does the 
result check for any special, simple value of (pi Ans: (voab tan (j)) / 3 

9.2 Repeat the preceding problem using the cylindrical surface of Eq. (9.4), but place the bottom point 
of the cylinder at coordinate (x, y) = (xq, 0). Ans: [v^a/ A){2xq + Tib/ A) 

9.3 Use the same velocity function v = xv^xy/b 2 and evaluate the flow integral outward from the 
closed surface of the rectangular box, (c < x < d), (0 < y < b), (0 < z < a). The convention is that 
the unit normal vector points outward from the six faces of the box. Ans: v^a{d — c)/2 

9.4 Work out the details, deriving the divergence of a vector field in spherical coordinates, Eq. (9.16). 

9.5 (a) For the vector field v = Af, that is pointing away from the origin with a magnitude proportional 
to the distance from the origin, express this in rectangular components and compute its divergence. 

(b) Repeat this in cylindrical coordinates (still pointing away from the origin though). 

(c) Repeat this in spherical coordinates, Eq. (9.16). 

9.6 Gauss's law for electromagnetism says § E-dA = q enc \/£o- If the electric field is given to be 

E = Ar, what is the surface integral of E over the whole closed surface of the cube that spans the 
region from the origin to (x,y,z) = (a, a, a)? 

(a) What is the charge enclosed in the cube? 

(b) Compute the volume integral, f d 3 rV -E inside the same cube? 

9.7 Evaluate the surface integral, §v-dA, of v = f Ar 2 sin 2 8 + 8 Br cos 8 sin <p over the surface of 
the sphere centered at the origin and of radius R. Recall section 8.8. 

9.8 (a) What is the area of the spherical cap on the surface of a sphere of radius R: < 8 < 8o? 

(b) Does the result have the correct behavior for both small and large 8q7 

(c) What are the surface integrals over this cap of the vector field V = fv o cos 8 sin 2 0? Consider both 
Jv-dA and / v x dA. Ans: v nR 2 (l - cos 2 8 )/2 

9.9 A rectangular area is specified parallel to the x-y plane at z = d and 0<x<a, a < y < b. A 
vector field is v = (xAxyz + yByx 2 + 2Cx 2 yz 2 ) Evaluate the two integrals over this surface 

J v-dA, and J dA x v 

9.10 For the vector field v = Ar n r, compute the integral over the surface of a sphere of radius R 
centered at the origin: <fv-dA. (n > 0) 

Compute the integral over the volume of this same sphere J d 3 rV -v. 

9.11 The velocity of a point in a rotating rigid body is v = tu x r. See problem 7.5. Compute its 
divergence and curl. Do this in rectangular, cylindrical, and spherical coordinates. 

9.12 Fill in the missing steps in the calculation of Eq. (9.29). 
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9.13 Mimic the calculation in section 9.6 for the divergence in cylindrical coordinates, computing the 
curl in cylindrical coordinates, V x v. Ans: Eq. (9.32). 

9.14 Another way to get to Eq. (9.39) is to work with Eq. (9.38) directly and to write the function 
p(r) explicitly as two cases: r < R and r > R. Multiply Eq. (9.38) by r 2 and integrate it from zero 
to r, being careful to handle the integral differently when the upper limit is < R and when it is > R. 



Note: This is not simply reproducing that calculation that I've already done. This is doing it a different 
way. 

9.15 If you have a very large (assume it's infinite) slab of mass of thickness d the gravitational field 
will be perpendicular to its plane. To be specific, say that there is a uniform mass density po between 
z = ±d/2 and that g = g z {z)z. Use Eqs. (9.36) to find g z {z). 

Be precise in your reasoning when you evaluate any constants. (What happens when you rotate the 
system about the x-axis?) Does your graph of the result make sense? 
Ans: in part, g z = +27iGpod, (z < —d/2) 

9.16 Use Eqs. (9.36) to find the gravitational field of a very long solid cylinder of uniform mass density 
po and radius R. (Assume it's infinitely long.) Start from the assumption that in cylindrical coordinates 
the field is g = g r {r, <fi, z)f, and apply both equations. 

Ans: in part g r = —2irGpor, (0 < r < R) 

9.17 The gravitational field in a spherical region r < R is stated to be g(r) = — fC/r, where C is a 
constant. What mass density does this imply? 

If there is no mass for r > R, what is g there? 

9.18 In Eq. (8.23) you have an approximate expression for the gravitational field of Earth, including 
the effect of the equatorial bulge. Does it satisfy Eqs. (9.36)? (r > i?Earth) 

9.19 Compute the divergence of the velocity function in problem 9.3 and integrate this divergence over 
the volume of the box specified there. Ans: {d — c)avo 

9.20 The gravitational potential, equation (9.42), for the case that the mass density is zero says to 
set the Laplacian Eq. (9.43) equal to zero. Assume a solution to V 2 V = to be a function of the 
spherical coordinates r and 9 alone and that 



Show that this works provided that / satisfies a certain differential equation and show that it is the 
Legendre equation of Eq. (4.18) and section 4.11. 

9.21 The volume energy density, dll / dV in the electric field is e^E 2 /2. The electrostatic field equations 
are the same as the gravitational field equations, Eq. (9.36). 




V(r,9) = Ar~( i+ Vf(x), 



where 



x = cos 6 



V-E = p/e 



and 



V x E = 



A uniformly charged ball of radius R has charge density po for r < i?, Q = 47rpo-R 3 /3. 
(a) What is the electric field everywhere due to this charge distribution? 
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(b) The total energy of this electric field is the integral over all space of the energy density. What is it? 

(c) If you want to account for the mass of the electron by saying that all this energy that you just 
computed is the electron's mass via Eq = mc 2 , then what must the electron's radius be? What is its 
numerical value? Ans: r e = 3 /§ (e 2 / Ane^inc 2 ) = 1.69fm 

9.22 The equations relating a magnetic field, B, to the current producing it are, for the stationary 
case, 

V x B = fj, J and V ■ B = 

Here J is the current density, current per area, defined so that across a tiny area dA the current that 
flows through the area is dl = J • dA. A cylindrical wire of radius R carries a total current / distributed 
uniformly across the cross section of the wire. Put the z-axis of a cylindrical coordinate system along 
the central axis of the wire with positive z in the direction of the current flow. Write the function J 
explicitly in these coordinates (for all values of r < R, r > R). Use the curl and divergence expressed in 
cylindrical coordinates and assume a solution in the form B = B \p(r , <fi , z) . Write out the divergence 

and curl equations and show that you can satisfy these equations relating J and B with such a form, 
solving for Ba. Sketch a graph of the result. At a certain point in the calculation you will have to 

match the boundary conditions at r = R. Recall that the tangential component of B (here Ba) is 

continuous at the boundary. 

Ans: in part, p Ir/2irR 2 (r < R) 

9.23 A long cylinder of radius R has a uniform charge density inside it, po ar, d it is rotating about its 
long axis with angular speed u. This provides an azimuthal current density J = porujcj) in cylindrical 
coordinates. Assume the form of the magnetic field that this creates has only a 2-component: B = 
B z (r,(J),z)z and apply the equations of the preceding problem to determine this field both inside and 
outside. The continuity condition at r = R is that the tangential component of B (here it is B z ) is 
continuous there. The divergence and the curl equations will (almost) determine the rest. Ans: in part, 
-pr 2 uj/2 + C (r < R) 

9.24 By analogy to Eqs. (9.9) and (9.17) the expression 

lim — 
V-*) V 

is the gradient of the scalar function 0. Compute this in rectangular coordinates by mimicking the 
derivation that led to Eq. (9.11) or (9.15), showing that it has the correct components. 

9.25 (a) A fluid of possibly non-uniform mass density is in equilibrium in a possibly non-uniform 
gravitational field. Pick a volume and write down the total force vector on the fluid in that volume; the 
things acting on it are gravity and the surrounding fluid. Take the limit as the volume shrinks to zero, 
and use the result of the preceding problem in order to get the equation for equilibrium. 

(b) Now apply the result to the special case of a uniform gravitational field and a constant mass density 
to find the pressure variation with height. Starting from an atmospheric pressure of 1.01 x 10 5 N/ m , 
how far must you go under water to reach double this pressure? 
Ans: Vp = —pg; about 10 meters 

9.26 The volume energy density, u = dll/dV, in the gravitational field is g 2 /8nG. [Check the units 
to see if it makes sense.] Use the results found in Eq. (9.39) for the gravitational field of a spherical 
mass and get the energy density. An extension of Newton's theory of gravity is that the source of 
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gravity is energy not just mass! This energy that you just computed from the gravitational field is then 
the source of more gravity, and this energy density contributes as a mass density p = u/c 2 would. 

(a) Find the additional gravitational field g r {f) that this provides and add it to the previous result for 

9t(t). 

(b) For our sun, its mass is 2 x 10 30 kg and its radius is 700,000 km. Assume its density is constant 
throughout so that you can apply the results of this problem. At the sun's surface, what is the ratio of 
this correction to the original value? 

(c) What radius would the sun have to be so that this correction is equal to the original g r (R), resulting 
in double gravity? Ans: (a) gcorrection/sWiginai = GM/lORc 2 

9.27 Continuing the ideas of the preceding problem, the energy density, u = dU/dV, in the gravita- 
tional field is g 2 /8nG, and the source of gravity is energy not just mass. In the region of space that is 
empty of matter, show that the divergence equation for gravity, (9.36), then becomes 

V -g = -4tcGu/c 2 = -g 2 /2c 2 

Assume that you have a spherically symmetric system, g = g r (r)f, and write the differential equation 
for g r . 

(a) Solve it and apply the boundary condition that as r -> oo, the gravitational field should go to 
g r (r) —GM/r 2 . How does this solution behave as r — > and compare its behavior to that of the 
usual gravitational field of a point mass. 

(b) Can you explain why the behavior is different? Note that in this problem it is the gravitational field 
itself that is the source of the gravitational field; mass as such isn't present. 

(c) A characteristic length appears in this calculation. Evaluate it for the sun. It is 1/4 the Schwarzchild 
radius that appears in the theory of general relativity. 

Ans: (a) g r = -GM/[r{r + R)] , where R = GM/2c 2 

9.28 In the preceding problem, what is the total energy in the gravitational field, f udV? How does 
this (-f-c 2 ) compare to the mass M that you used in setting the value of g T as r — > oo? 

9.29 Verify that the solution Eq. (9.48) does satisfy the continuity conditions on V and V. 

9.30 The r-derivatives in Eq. (9.43), spherical coordinates, can be written in a different and more 
convenient form. Show that they are equivalent to 

1 d 2 (rV) 
r Or 2 

9.31 The gravitational potential from a point mass M is —GM/r where r is the distance to the mass. 
Place a single point mass at coordinates (x,y,z) = (0, 0, d) and write its potential V. Write this 
expression in terms of spherical coordinates about the origin, (r, 6), and then expand it for the case 
r > d in a power series in d/r, putting together the like powers of d/r. Do this through order (d/r) 3 . 
Express the result in the language of Eq. (4.61). 

Ans: -GU - GMd[ COS 8] - £Md! [ 3 COS 2 # - ±] _ GM^[5 cog 3 # _ 3 cog 0] 

9.32 As in the preceding problem a point mass M has potential —GM/r where r is the distance to 
the mass. The mass is at coordinates (x,y,z) = (0,0, d). Write its potential V in terms of spherical 
coordinates about the origin, (r, 9), but this time take r < d and expand it in a power series in r/d. 
Do this through order (r/d) 3 . 

Ans: {-GM/d)\l + (r/d)Pi(cos0) + (r 2 /d 2 )P 2 (cos#) + (r 3 /d 3 )P 3 {cos6) + • • •] 
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9.33 Theorem: Given that a vector field satisfies V x v = everywhere, then it follows that you can 
write v as the gradient of a scalar function, v = — Vijj. For each of the following vector fields find if 
possible, and probably by trail and error if so, a function ijj that does this. First determine is the curl is 
zero, because if it isn't then your hunt for a ip will be futile. You should try however — the hunt will 
be instructive. 

(a) xy 3 + 3yxy 2 , (c) xycos(xy) + yxcos(xy), 

(b) xx 2 y + yxy 2 , (d) xy 2 s'mh(2xy 2 ) + 2y xy smh(2xy 2 ) 

9.34 A hollow sphere has inner radius a, outer radius b, and mass M, with uniform mass density in 
this region. 

(a) Find (and sketch) its gravitational field g r {r) everywhere. 

(b) What happens in the limit that a — > bl In this limiting case, graph g r . Use g r {r) = —dV/dr and 
compute and graph the potential function V(r) for this limiting case. This violates Eq. (9.47). Why? 

(c) Compute the area mass density, a = dM/dA, in this limiting case and find the relationship between 
the discontinuity in dV/dr and the value of a. 

9.35 Evaluate 

^ij^ijki ^mjk^njki 9{Xj, 9{Xi, ^ijk^ijki ^ij^j 
and show that ^ijk e mnk = ^im^jn ~ ^in^jm 

9.36 Verify the identities for arbitrary A, 

(A-V)f=A or AidiXj = Aj 
V ■ V x v = or dieijkdjVk = 
V ■ (fA) = (Vf) -A + f(V-A) or d^fAA = (d^A, + fd^ 

You can try proving all these in the standard vector notation, but use the index notation instead. It's 
a lot easier. 

9.37 Use index notation to prove V x V x v = V(V ■ v ) — V 2 v. First, what identity you have to prove 
about e's. 

9.38 IsVxjJ perpendicular to V? Either prove it's true or give an explicit example for which it is 
false. 

9.39 If for arbitrary A{ and arbitrary Bj it is known that aijA^Bj = 0, prove then that all the aij are 
zero. 

9.40 Compute the divergences of 

Axx + By 2 y + C z in rectangular coordinates. 
Arf + B6 2 6 + C (f) in spherical coordinates. 
How do the pictures of these vector fields correspond to the results of these calculations? 

9.41 Compute the divergence and the curl of 

yx — xy , . yx — xy 

and of 



x 2 + y 2 ' (x 2 + y 2 ) 2 
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9.42 Translate the preceding vector fields into polar coordinates, then take their divergence and curl. 
And of course draw some pictures. 

9.43 As a review of ordinary vector algebra, and perhaps some practice in using index notation, translate 
the triple scalar product into index notation and prove first that it is invariant under cyclic permutations 
of the vectors. 

{a) A-B xC = B-C xA = C-Ax B. Then that 

(b) A-Bx C = AxB- C. 

(c) What is the result of interchanging any pair of the vectors in the product? 

(d) Show why the geometric interpretation of this product is as the volume of a parallelepiped. 

9.44 What is the total flux, § E • dA, out of the cube of side a with one corner at the origin? 

(a) E = ax + (3y + 7 z 

(b) E = axx + (3yy + •yzz. 

9.45 The electric potential from a single point charge q is kq/r. Two charges are on the z-axis: — q 
at position z = Zq and +q at position Zq + a. 

(a) Write the total potential at the point (r,9,(j)) in spherical coordinates. 

(b) Assume that r S> a and r 3> Zq, and use the binomial expansion to find the series expansion for 
the total potential out to terms of order l/r 3 . 

(c) how does the coefficient of the l/r 2 term depend on ZqI The coefficient of the l/r 3 term? These 
tell you the total electric dipole moment and the total quadrupole moment. 

(d) What is the curl of the gradient of each of these two terms? 
The polynomials of section 4.11 will appear here, with argument cos#. 

9.46 For two point charges q\ and qi, the electric field very far away will look like 
that of a single point charge q\ + 52- Go the next step beyond this and show that 
the electric field at large distances will approach a direction such that it points along 
a line that passes through the "center of charge" (like the center of mass): {q\f\ + 
<?2^*2)/(<?i + 92)- What happens to this calculation if q2 = — q{? You will find the 
results of problem 9.31 useful. Sketch various cases of course. At a certain point in the 
calculation, you will probably want to pick a particular coordinate system and to place 
the charges conveniently, probably one at the origin and the other on the z-axis. You 
should keep terms in the expansion for the potential up through l/r 2 and then take — VV. Unless of 
course you find a better way. 

9.47 Fill in the missing steps in deriving Eq. (9.58). 

9.48 Analyze the behavior of Eq. (9.58). The first thing you will have to do is to derive the behavior 
of sinh -1 in various domains and maybe to do some power series expansions. In every case seek an 
explanation of why the result comes out as it does. 

(a) If z = and r = ^/ x 2 + y 2 3> L what is it and what should it be? (And no, zero won't do.) 

(b) If z = and r <<C.L what is it and what should it be? 

(c) If z > L and r — > what is this and what should it be? Be careful with your square roots here. 

(d) What is the result of (c) for z > L and for z - L < LI 

9.49 Use the magnetic field equations as in problem 9.22, and assume a current density that is purely 
azimuthal and dependent on r alone. That is, in cylindrical coordinates it is 



J = Mf(r) 
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Look for a solution for the magnetic field in the form B = 2B z (r). What is the total current per 
length in this case? That is, for a length Az how much current is going around the axis and how is 
this related to B z along r = 0? 

Examine also the special case for which f(r) = except in a narrow range a < r < b with b — a <^ib 
(thin). Compare this result to what you find in introductory texts about the solenoid. 

9.50 For a spherical mass distribution such as the Earth, what would the mass density function have 
to be so that the gravitational field has constant strength as a function of depth? Ans: p oc l/r. 

9.51 Use index notation to derive 



9.53 Use index notation to compute Ve . Also compute the Laplacian of the same exponential, 
V 2 = div grad. 

9.54 Derive the force of one charged ring on another, as shown in equation (2.35). 

9.55 A point mass m is placed at a distance d > R from the center of a spherical shell of radius R 
and mass M. Starting from Newton's gravitational law for point masses, derive the force on m from 
M. Place m on the z-axis and use spherical coordinates to express the piece of dM within d6 and 
d<p. (This problem slowed Newton down when he first tried to solve it, as he had to stop and invent 
integral calculus first.) 

9.56 The volume charge density is measured near the surface of a (not quite perfect) conductor to 
be p(x) = p§e x l a for x < 0. The conductor occupies the region x < 0, so the electric field in the 
conductor is zero once you're past the thin surface charge density. Find the electric field everywhere 
and graph it. Assume that everything is independent of y and z (and t). 



(AxB)-(CxD) 



(A-C)(B-D) - (A-D)(B-C) 



9.52 Show that V -(A x B) = B -V x A 



A-V x B. Use index notation to derive this. 



